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Abstract 

We present a magnetometer based on optically pumped Cs atoms that 
measures the magnitude and direction of a 1 pT magnetic field. Multiple 
circularly polarized laser beams were used to probe the free spin preces¬ 
sion of the Cs atoms. The design was optimized for long-time stability 
and achieves a scalar resolution better than 300 fT for integration times 
ranging from 80 ms to 1000 s. The best scalar resolution of less than 80 fT 
was reached with integration times of 1.6 to 6 s. We were able to mea¬ 
sure the magnetic field direction with a resolution better than 10 prad for 
integration times from 10 s up to 2000 s. 
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1 Introduction 


Magnetometers using optical pumping (OPM) of atomic media were pioneered 
in the early 1960s [T]. Since then, many OPM varieties [2] have been developed 
for diverse applications, e.g., mapping the geo-magnetic field or detection of the 
bio-magnetic field emanating from the human heart EE and brain EE]. In 
fundamental science, OPMs monitor the magnetic field in precision magnetic 
resonance experiments searching for electric dipole moments (EDM) E E and 
Lorentz invariance tests EES- The neutron EDM (nEDM) search sets strin¬ 
gent constraints on theories proposing extensions beyond the standard model of 
particle physics. Experimental sensitivity to a nEDM depends directly on the 
control and measurement of the magnetic field in the experiment, a task of par¬ 
ticular challenge since (currently) the field must be known over volumes on the 
order of 20 1 for times of hundreds of seconds. Herein, we present an OPM com¬ 
bining long-term stability with high statistical sensitivity, and including vector 
information. The OPM is designed to serve in an array of such sensors to form 
an auxiliary magnetometer system monitoring the stability and uniformity of 
the magnetic field in a next generation nEDM experiment at the Paul Scherrer 
Institut. An array of scalar Cs OPMs has been used successfully to determine 
directional and gradient magnetic field information via fits to multi-sensor read¬ 
ings fm . 

Various methods exist for extracting information about the magnetic field 
vector components, including the spin exchange relaxation free (SERF) mag¬ 
netometers E operating at | Bo | = 0 and whose intrinsically sensitivity is to 
one vector component only. For operation in the offset fields used for neutron 
magnetic resonance, we focus on conventional OPMs that measure the magnetic 
field modulus by detecting the Larmor precession frequency, wl = 7 |B 0 ], where 
7 is the gyromagnetic ratio of the probed atomic state. Information about vec¬ 
tor components is gleaned by monitoring the OPM’s response to an externally 
applied oscillating magnetic field using phase sensitive detection. In a recently 
published all-optical variant of that method [12] , circularly polarized laser beams 
induce an effect equivalent to the perturbation field via the vector light shift 
jl'3] . Without external modulations, vector information can be inferred using 
multiple detection channels m when the first and second harmonics of the Lar¬ 
mor precession of an atomic alignment m are detected with linearly polarized 
light. 

Our approach uses multiple circularly polarized laser beams to gain vector 
information and thus extends methods pioneered by Fairweather and Usher DU- 
The absorption of circularly polarized light depends linearly on the projection 
of the atomic spin polarization on the light’s k vector. The precessing atomic 
polarization modulates the transmitted light power at ui if Bo is not parallel 
to k. This can be used to maintain the condition Bo || k in a feedback loop 
either by changing the direction of Bo [T6] or of k m■ In contrast to those 
vector magnetometer implementations, our system uses off-line data analysis 
enabling us to infer the magnetic field information from free spin precession 
(FSP) signals. The FSP method is particularly well suited for the application 
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Figure 1: a) Apparatus schematic: four laser light beams are polarized by 
the optical modules (om), traverse the Cs cell (ce, shown in section) and are 
converted into signals 5) by the photodiodes (pd). b) Measurement cycle time 
structure, repeated every 40 ms: The longitudinal magnetization M\ is created 
by optical pumping. Following a 7t/2 pulse, B y , of duration qT 2 , the resultant 
transverse magnetization M t gives rise to FSP signals (shown here projected 
onto the x axis). The parameters a, /3 measure the length of the n/2 pulse and 
the FSP signal, respectively, in units of T 2 . In each signal analysis, the t = 0 
time origin is reset to the FSP start. 


in nEDM experiments since it allows for very stable field measurements. 


2 Experimental setup 

The experiment was performed inside the magnetic shield of the nEDM experi¬ 
ment at PSI [18J, in which a stable 1 pT magnetic field was generated by a cos 9 
coil. The static magnetic field is parametrized as Bo = I?o(sm 9 cos <p, sin 9 sin </>, 
cos 9), and was approximately aligned along the z axis, i.e., 9 r; 0. Our mag¬ 
netometer was created to measure the field modulus, Bq, and its direction, i.e., 
the polar angle, 9 , and azimuthal angle <j). The magnetometer design is shown 
in Fig. Oja). Light is generated by an extended cavity diode-laser coupled to a 
polarization-maintaining single-mode fiber splitter having three outputs. One 
output feeds a saturated absorption spectroscopy unit, used for active laser fre¬ 
quency stabilization to the F =4 —> 3 cesium D\ transition (894 nm). A second 
single-mode fiber guides light to the magnetometer head, where it is split into 
four beams which are coupled into short multi-mode fibers. At the sensor head, 
the light from each multi-mode fiber is collimated and circularly polarized by a 
linear polarizer and a quarter-wave plate, mounted in a compact optical module 
(om). The power of those beams can be adjusted by rotating an additional linear 
polarizer in the om. The beams traverse an evacuated 45 mm diameter glass-cell 
(ce) containing a saturated vapor of cesium atoms. The cell is paraffin coated 
[l9l to reduce spin depolarization during atom-wall collisions. A combination 
of photodiodes (pd) and transimpedance amplifiers converts the transmitted 
light power of each beam to a signal S), which is digitized with a high resolution 
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sampling system. The combined noise of the photodiode, preamp, and sampling 
system is well below the shot-noise level for the typical light power of 1 pW per 
laser beam. 

The magnetometer is operated in pulsed mode, and information is extracted 
from the FSP signals. Figure [IJb) shows the experimental cycle which repeats 
every 40 ms. The FSP is described using the magnetization M associated with 
the ensemble average of the atomic spin. The combined optical pumping by 
the four laser beams in combination with Bo creates a magnetization Mi lon¬ 
gitudinal to Bo- A short magnetic 7r/2 pulse along the y direction turns Mi to 
a direction transverse to Bo, creating M t . The pulse uses a single sinusoidal 
period in order to minimize deadtime (a in Fig. [T]b)). The magnetization 
component perpendicular to Bo precesses at the Larmor frequency, ul = 7 B 0 , 
where 7 = 3.498 6211(4) kHz/pT is the gyromagnetic ratio of the F =4 cesium 
ground state. The light absorption by the cesium atoms depends linearly on 
the projection of M on the light’s k -vector j2U]. Consequently, the transmitted 
laser power, measured by a photodiode, is modulated at u l- The transverse 
magnetization component decay (see M t in Fig. [lj)b)), with its effective decay 
time T 2 , is observed as a decreasing modulation amplitude of the recorded FSP 
photodiode signal. During the FSP, the longitudinal magnetization, M\, is recre¬ 
ated by optical pumping, such that the next n/2 pulse can start the next cycle. 
Both the data acquisition system recording the FSP signals and the function 
generator producing the n/2 pulses are synchronized to an atomic clock. 

Using the classical Bloch equation, the recorded signal for each laser beam 
can be modeled as 


_ t 

Si(t)=Ci + e t 2 ( 6 j —i— A Ci i cosut + A Sii sinwt) . (1) 

Both frequency ui and effective decay time T 2 are common parameters for all 
simultaneously recorded FSP signals. The offsets, c t and bi , as well as the in- 
phase, A c> i, and quadrature, A s> j, components of the modulation amplitudes 
are different for each signal Si- The Ci parameters represent the DC signal 
offsets and are proportional to the average light power of beam i. If the k- 
vector of beam i has a longitudinal component, the exponential build-up of 
M\ contributes to the absorption it probes. Assuming that the longitudinal and 
transverse relaxation rates are equal allows this contribution to be parametrized 
by the offsets bi. The modulation amplitudes are used to determine the magnetic 
field direction. The Bo field magnitude is determined using the estimation of 
frequency w, interpreted as the Larmor frequency. The magnitude and the 
extracted field direction are used to reconstruct the vector magnetic field. 

3 Data analysis 

The parameters of Eq. © are extracted with a precision limited by the Cramer- 
Rao Lower Bound (CRLB) [21] [22.. The lower limit of the frequency spectral 
density, py , calculated with the CRLB for signals with no DC components (c, = 
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Figure 2: The Larmor precession is described by M c and M s , c.f. Eq. J3]), 
which together define a plane perpendicular to Bo (here depicted as a circle. 
The longitudinal magnetization Mi changes only slowly due to relaxation and 
optical pumping. Signal Si(t) represents the transmitted power of the laser 
beam k;. 


bi = 0) and sampled at a sufficiently high rate (3> wl/27t) is 

2 p I (a + /3)eP sinh/3 

Pf ~ ttAT 2 \I cosh 2/3-2/3 2 -l ' U 

The length of the FSP signal, T, is parametrized in a dimensionless way by 
/? = T/Tn (Fig. [Ub)), whereas a measures the dead-time of the 7r/2 pulse. The 
spectral density p of the photodiode signals is ultimately limited by shot noise. 
The amplitude A is proportional to M z , in the instance before it is flipped, 
it thus scales like A = Aq (l — e~ ^). Given this, Eq. ([2]) has a minimum at 
/? ~ 2.6. For technical reasons, the pulse repetition time of T = 40 ms was 
chosen to be slightly shorter than the optimum given that = 20.4(2) ms. 
Estimates indicate a ~ 10% performance gain at Tk 2.6 T\ 

Reconstructing the Bo vector components is made by monitoring the M 
component precessing at wt. By definition, the precession happens in a plane 
perpendicular to Bo, thus the cross product of two vectors in that plane yields 
a vector parallel to Bo- The method’s statistical sensitivity is maximized when 
the phase difference of the two vectors is 7r/2. This is achieved by parametrizing 
the precessing part of M by its in-phase and quadrature components as 

Mt(i) = Me cos wi A + M s sin wiT • (3) 

If the two vectors are known, the Bo direction follows as 

B 0 oc M c x M s . (4) 

Measuring M c and M s is straightforward in a three-beam magnetometer with 
orientations along the Cartesian coordinates axes (Fig. [2]). In such a configu¬ 
ration, the cos and sin modulation components seen by each beam (Eq. Jl])) 
correspond directly to M c and M s . 
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The experimental vector magnetometer reported herein uses four laser beams 
(Fig. Ufa)). In this configuration, each beam has a k component along z , thus 
contributing to optical pumping provided Bo is approximately oriented along 
z, which is our case. Each beam probes the projection, Mj = fcj • M, of the 
magnetization M onto its /c-vector. The 3D vector M is reconstructed from the 
four projections, using the projection matrix P 



( Ml \ 

M 2 

m 3 

\m 4 ) 





(5) 


This reconstruction is advantageous since M x and M y are derived by subtract¬ 
ing two projections, reducing common mode noise. Since M x and M y are the 
important components for determining the Larmor frequency, the chosen beam 
configuration facilitates its high resolution extraction. All parts of the signal 
Si(t) (Eq. (JTJ)) that depend linearly on M are transformed from the four projec¬ 
tions into a 3D representation using the matrix P. In the low light-power limit, 
the in-phase and quadrature modulation amplitudes of 5,(t) are proportional to 
the DC signal, Cj, detected by photodiode i. Amplitudes A c i and A s i are nor¬ 
malized using Cj, to compensate for slight differences in light power and possible 
differences in photodiode preamplification factors. Finally, P converts the nor¬ 
malized amplitudes extracted from signals S/(f) to 3D vectors that determine 
the direction of B 0 using Eq. ©. Using numerical simulations |23l . we veri¬ 
fied that the resulting angles 9 and <j> are determined with maximum statistical 
efficiency. The magnetic field direction and magnitude can be extracted from 
the data of only three laser beams. Given the geometry of the beams in the 
experiment this is, however, not possible at maximum statistical efficiency. Cor¬ 
relations in the signals S t (t) due to the over determined measurement with four 
laser beams can be used to verify the normalization factors of the amplitudes 


The projection matrix P depends on the actual orientation of the laser 
beams. Deviations from the assumed orientations lead to systematic errors 
in the extracted magnetic field orientation 9 and 4>. Those errors depend in a 
complex way on the orientation of the magnetic field and the direction in which 
the beam is tilted. If one laser beam is tilted by an angle A a in a direction that 
causes the largest errors, it contributes an error of A 9 = 1/4 Aa to the extracted 
magnetic field orientation. Tilting all four laser beams in this way is equivalent 
to tilting the whole sensor by Aa which naturally causes an estimation error 
of A 9 = A a. Tilting all beams in random directions causes a combined error 
of A 9 = 1/2 Aa. The mechanical construction of the experiment can currently 
not guarantee a alignment better than Aa = 0.004 rad. 

Two estimation methods to extract the parameters of Eq. m from the digi¬ 
tized signals were studied: Least-squares fitting, and demodulation. The least- 
squares method fits the Eq. m model to the experimental data gained from all 
beams simultaneously. 
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Figure 3: ASD of the magnetic field magnitude (a, b), the field’s z-component 
(c, d), and the B direction 6 (e, f) measured with the vector cesium magnetome¬ 
ter. The results obtained by least-squares fitting (curves a, c, and e) generally 
show better ASD values at short integration times, while the demodulation 
method (curves b, d, f) provides better results over longer integration times. 
Curves c and e are affected by systematic errors for integration times larger 
than 4 s. The error bars were calculated according to (25] . 


The demodulation method uses two-phase lock-in detection with cos and 
sin reference signals at a frequency w r close to wl- This mixes the w l mod¬ 
ulation down to a frequency close to DC, while noise and other modulations 
are suppressed by the low-pass filter [24] . The in-phase and quadrature lock-in 
signals are converted to phase <p(tj) = ipj and amplitude A(tj) = Aj for each 
FSP signal. The initial modulation amplitude, A(t = 0), is extracted by a 
least-squares fit of A(0) exp —t/T 2 to the time series Aj of one FSP. The model 
ip(t) = ip(0)+ojt is fitted (with weighting factors 1 /Aj) to the phase signal after 
correcting for discrete 2n steps. The frequency difference between w r and wl is 
found using the slope u>. 

The in-phase and quadrature modulation amplitudes are found via A C j = 
Acostp(O), and A S) j = — Asin<p(0). For both methods, the least-squares fits are 
made simultaneously for the four FSP signals using one common frequency, w, 
and decay time, T?. 

4 Results 

The nEDM experiment requires magnetic field measurements that are stable 
over hundreds of seconds at the sub-pT level. The cesium vector magnetome- 
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ter statistical errors, as described by the CRLB, are by far sufficient to reach 
that goal. However, systematic errors arising from drifting parameters limit the 
long-term stability that this magnetometer can achieve. To characterize the 
long-term stability, we measured during 11 hours under best-case conditions of 
field stability. Figure [3] shows the Allan standard deviation (ASD) [5Bj of vector 
and scalar field measurements as a function of integration time r. Using the 
estimated FSP parameters and a noise density extracted from the measured 
data’s Fourier spectrum, Eq. Pjl yields a CRLB of 81 fT/ -\/Hz for the held 
magnitude. For r < 1 s, the ASD plots show the expected improvement propor¬ 
tional to t -1 / 2 . The least-squares fitting possesses a higher statistical efficiency, 
which is visible as a 12% smaller ASD. This difference disappears for r > Is 
where the ASD is no longer limited by statistical processes. For the longer in¬ 
tegration times, the ASD is limited by magnetic held drifts and magnetometer 
instabilities, thus, the ASD represents the limit for the magnetometer stability. 
The magnitude can be measured with an uncertainty smaller than 300 fT for 
integration times ranging from 80 ms to 1000 s. A best sensitivity of 75.2 fT is 
achieved at r = 4s, which corresponds to a relative sensitivity of 7.6 x 10 -8 . 

Curves c and d in Fig.[3]show the ASD for the held’s z component (longitudi¬ 
nal) measurement. This shows that the values extracted using the demodulation 
method are more stable than those from least-squares htting for long integration 
times. This happens because the least-squares htting does not model the second 
harmonic of the Larmor modulation, 2 wl, while the demodulation method is 
insensitive to it due to the low-pass hlter. The integration times for which the 
z component can be measured with an uncertainty smaller than 300 fT range 
up to 1000 s, but start at 2 s, due to the larger statistical errors. This increase 
in statistical uncertainty is due to using amplitudes which cannot be estimated 
as precisely as the Larmor frequency. 

Figure [3] e and f show the ASD of 9, characterizing the direction Bq, as 
derived from the estimated vector components. The ASD of 9 , estimated using 
least-squares htting, scales statistically for r < 1 s and complies very well with 
the CRLB calculated using error propagation from the estimated amplitudes. 
Using demodulation, the resolution of 9 reaches 6.4prad for t = 10 s, and does 
not change significantly until r = 2000 s. The ASD of 4> behaves similarly, but 
with larger uncertainties since the measurement was made near the degenerate 
case of 9=0. 


5 Conclusion and discussion 

The presented magnetometer achieves high sensitivity both in magnitude and 
held direction. Upper limits on processes that limit the stability of the magne¬ 
tometer readings are derived from the ASD plots and show high sensitivity is 
maintained even at integration times of 1000 s. This value is probably limited 
by drifts of the Bq held components in the present nEDM experiment. Further 
studies will try to distinguish between instabilities intrinsic to the magnetometer 
and external held drifts by using several magnetometer modules. 
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In contrast to other recently published vector magnetometers m the pre¬ 
sented approach does not degrade the scalar resolution when extracting vector 
information. Consequently, it achieves an order of magnitude better scalar 
resolution while being able to resolve the Bg direction with 59 < lOprad for 
integration times ranging beyond 2000 s. This makes the presented approach 
an ideal choice for applications that use long integration times. For our future 
nEDM apparatus it is planned to use an array of vector Cs magnetometers in 
order to monitor the Bq field and its gradients. Scaling to multiple sensors is 
aided by the low needs on laser power and the efficient data processing possible 
in the demodulation mode. 

The magnetometer presented here requires calibration in order to provide 
absolute field direction information. However, the accuracy of its absolute field 
magnitude information may be limited—as discussed by Grujic et al in Eil¬ 
at the several 10 pT level since k JL Bq. The demonstrated stability at long 
integration times is a necessary step for the future development of such cali¬ 
bration procedures. With the stability proven, the detailed studies of device 
construction systematics (e.g., perturbations to the values in P (Eq. ©) and 
device alignment to an external coordinate system will permit calibration of the 
device, thus moving it from being a field stability measurement system to a true 
field measurement system. 

A remaining disadvantage of this approach is the non-‘magnetically silent’ 
7 t/ 2 spin manipulation pulse, which can perturb the environment under study. 
A straightforward way to overcome this is the use of Bell-Bloom pumping, cur¬ 
rently under development within our collaboration |27] . A combination of these 
two methods is being pursued to provide a sensitive and magnetically silent 
vector magnetometer for our future nEDM search. 
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